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"q^" ■ It has been proposed that an asymmetry in the electron neutrino sector 

On 

Q^ ■ may be generated by resonant active-sterile neutrino transformations during 

£>y Big Bang Nucleosynthesis (BBN). We calculate the change in the primordial 

a 

4 He yield Y resulting from this asymmetry, taking into account both the time 



(N 



evolution of the v e and v e distribution function and the spectral distortions 



fN) . in these. We calculate this change in two schemes: (1) a lepton asymmetry 

. directly generated by v e mixing with a lighter right-handed sterile neutrino 

. v s \ and (2) a lepton asymmetry generated by a v T «-» v s or «-> u s trans- 

o. 

s * ■ formation which is subsequently partially converted to an asymmetry in the 

. z/ e i/ e sector by a matter-enhanced active-active neutrino transformation. In 

q . the first scheme, we find that the percentage change in Y is between —1% 

and 9% (with the sign depending on the sign of the asymmetry), bounded 
by the Majorana mass limit m Ve < 1 eV. In the second scheme, the max- 
imal percentage reduction in Y is 2%, if the lepton number asymmetry in 
neutrinos is positive; Otherwise, the percentage increase in Y is < 5% for 



o 



m?, „ — m 



1 < 10 4 eV. We conclude that the change in the primordial 4 He 



yield induced by a neutrino-mixing-generated lepton number asymmetry can 
be substantial in the upward direction, but limited in the downward direction. 

PACS numbers: 14.60.Pq; 14.60.St; 26.35. +c; 95.30.-k 
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I. INTRODUCTION 



It has been known that a resonant active-sterile neutrino transformation during the Big 
Bang Nucleosynthesis (BBN) epoch can generate lepton number asymmetries in the active 
neutrino sectors [ED-fl. The generated lepton number asymmetry L Va (with v a being any of 
the three active neutrino species) has an order of magnitude 

t __ n v a ~ n v a , \Sm 2 /eV 2 \ 



for \L Va \ < 0.1. (1) 



n 7 10 (T/MeV) 4 

In this equation, n is the particle proper number density, 5m 2 = m? Us — m? v , and T is the 
temperature of the universe. 

The asymmetry may have an appreciable impact on the the primordial 4 He abundance Y 
if it is in the v e sector and if its magnitude at the weak freeze-out temperature T ~ 1 MeV 
is |Lj,J > 0.01. The resulting change of the primordial 4 He abundance, AY, however, is not 
easy to estimate. Not only is the generated asymmetry a function of time, but also the u e or 
v e energy spectrum is distorted by mixing (in a time- dependent fashion as well). Therefore, 
previous attempts || to estimate AY in this lepton number generation scenario, employing 
BBN calculations based on a constant asymmetry and a thermal neutrino spectrum for 
electron-type neutrinos, is overly simplistic and may yield inaccurate results. 

In this paper, we discuss in detail the time evolution of L Ue and the distortion in the v e or 
P e spectrum in this neutrino-mixing-driven lepton number asymmetry generating scenario. 
We then calculate AY taking into account these time-dependent and energy- dependent 
effects by modifying the standard BBN code accordingly. We consider two schemes: a direct 
one and an indirect one. 

The direct scheme involves a direct resonant v e <-> v s transformation which generates 
L Ve mil. The indirect one first has a lepton asymmetry generated via a resonant <-> v s 
or v T <-> v s transformation and then has the asymmetry partially transferred into L Ve by an 
active-active <-> v e or v T <-> v e transformation ||. In both cases, Eq. (P indicates that the 
active-sterile channel requires m 2 u —m 2 u > 0.1 eV 2 to generate \L Va \ > 0.01 at a temperature 
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T ~ 1 MeV. This implies m Va > 0.3 eV. In addition, the effective mixing angles associated 
with these neutrino oscillation channels have to be large enough to generate lepton number 
asymmetries efficiently, but not so large as to produce too many z/ s 's before or during the 
onset of the process. The excess z/ s 's may not only suppress the lepton number generation, 
but also yield a 4 He mass fraction that is too large These conditions can be quantified 
as: 

|5m 2 /eV 2 | 1/6 sin 2 2^> KT 11 , (2) 

\5m 2 /eV 2 \ sin 4 26 < 1(T 9 (1(T 7 ) for v a = u e v T ). (3) 

where 9 is the vacuum mixing angle. 

The change in the predicted primordial 4 He abundance can go either way, depending 
on whether L Ve is positive (decreasing Y) or negative (increasing Y). Our detailed BBN 
calculations show that in the direct scheme, we have —0.002 < AY < 0.022, bounded by the 
v e Majorana mass limit m Ue < 1 eV. This possible change in Y is significant compared to 
the uncertainty involved in the Y measurements. (The measured primordial 4 He abundance 
from Olive et al. is Y = 0.234 ±0.002 (stat.) ±0.005 (syst.) [[?]], while another group claims 
Y = 0.244 ±0.002 ||.) In fact, an increase in Y of magnitude > 0.01 due to a large negative 
L Ve would already be inconsistent with observations. In the indirect scheme, we find that 
the maximal possible reduction in Y is 0.005. The expected increase of Y goes up with 
m?, ,j — m?, > and can be as high as 0.013 for m 2 ,. — m 2 ~ 10 4 eV 2 . In both schemes, 
AY is rather limited in the negative direction. The possible reduction in Y, resulting from 
a positive L Ue , may to some degree narrow the gap between the lower Y measurement @ 
and the standard BBN prediction (Y = 0.246 ± 0.001 @ assuming a primordial deuterium 
abundance D/H^ 3.4 ± 0.3 x 10~ 5 ||10|| ). We note, however, that the maximal reduction of 
AF pa —0.005 is achieved only in the indirect scheme when m 2 .. — m 2 ~ 100 to 300 eV 2 , 
implying an unstable or v T with m v Vr > 15 eV 0. 
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II. GENERATION OF LEPTON ASYMMETRY BY RESONANT 
ACTIVE-STERILE NEUTRINO TRANSFORMATION 



The formalism of active-sterile neutrino transformation and associated amplification of 
lepton asymmetry when 5m 2 = m? v — m 2 a < has been discussed extensively elsewhere 
[ if-pUS, n|. Here we summarize the main conclusions of these papers, and then concentrate 



our discussion on the final stage of the amplification process, when T approaches the weak 
freeze-out temperature ~1 MeV. 

For v a <-> v s transformation with 5m 2 < 0, a resonance occurs at 

1/6 

(4) 



^res ~ Tq , 



^ -1/3 


5m 2 cos 28 




leV 2 



where T ~ 19(22) MeV for a = e (/i, r), and E/T is the neutrino energy normalized by the 
ambient temperature. Below T res , lepton asymmetry may be amplified to asymptotically 
approach one of the two values in eq. ([!]). Before this asymptotic value is reached, however, 
there is a brief "chaotic" phase in which L Va oscillates around zero [[|. As a result, the 
sign of L Ua that emerges from the chaotic phase is unpredictable. We should point out 
that the detailed numerical evolution of the generated lepton number remains controversial. 
For example, whether or not the evolution of the lepton number represents true chaos is 
not precisely known. However, our BBN arguments simply reply on the sensitivity of the 
generated lepton number to the neutrino oscillation parameters. This sensitivity leads to 
a causality consideration on the sign of the generated lepton number. A discussion of the 
causality consideration can be found in Ref. [12 . 



It is not surprising that at these two asymptotic L Va values, either v a <-> v s (if L Ua < 0) 
or v a <-> v s (if L Va > 0) undergoes resonant transition due to matter effects. The system 
maintains the growth of L Va by converting one of u a /i> a resonantly but suppressing the 
transformation of the other. Of course, not all v a or u a undergo the resonant transforma- 
tion, because neutrinos in the early universe have an energy distribution and the resonance 
condition is energy dependent. When L Ua is small, even the resonant conversion of a small 
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fraction of either u a or v a with energy E TeB will be enough to maintain the growth of L Ua . 
To quantify the above arguments, we note that the effective potential V= (V x , V y , V z ) of 
the v a <-> u s transformation channel is 

V x = -—sm28, ^ = 0, V^ = -— cos2fl + ^ + V; T . (5) 



The contribution from matter-antimatter asymmetries (matter effect) is [pl| 



V^±0.35G F T 3 [L + 2L„ a + £ L Vp ], (6) 

where Gp is the Fermi constant, and L ~ 10~ 9 represents the contributions from the 
baryonic asymmetry as well as the asymmetry in electron-positions. The "+" sign is for the 
neutrino oscillation channel, and the "— " sign is for the anti-neutrino oscillation channel. 
The contribution from the thermal neutrino background is Vj, whose value is 



v t w _ a V^l±J^ gI ET \ ( 7 ) 
where A w 110(30) for a = e (// or r). The effective matter mixing angle at temperature T 



is 



Kb 

sm2fl eff = (V r g2 + ^2)i/2 » ( 8 ) 

which reduces to vacuum mixing when l/f' and Vj^ are zero. 

Several physical processes with different times scales are involved in the resonant v a <-> z/ s 
transformation process: (1) the local neutrino oscillation rate |V|; (2) the weak interaction 
rate T ~ 4G^T 5 ; (3) the Hubble expansion rate H = —T/T = 5.5T 2 /m p \ where m p \ ~ 
1.22 x 10 28 eV is the Planck mass; (4) the rate of change of |V|, |V|/|V|, caused by the 
change of lepton asymmetry and the Hubble expansion. If any one of the rates is much larger 
than the others, we may consider all the other processes as perturbations, which simplifies 
the picture greatly. For example, if | V| dominates we may consider the system as an ordinary 
neutrino transformation system with an effective mixing angle as in Eq. and with all 
physical variables changing adiabatically. If the weak interaction dominates, each weak 



scattering acts as a "measurement" to the mixing system, effectively reducing the mixture 
to flavor eigenstates u a and v s . The neutrino transformation is hence suppressed, with a 
reduced sin2# e g- = V x /(T/2). If the Hubble expansion dominates, the other processes are 
essentially "frozen out" . This is the case at T < 1 MeV when the two-body weak interaction 
freezes out and the neutron-to-proton ratio becomes fixed (other than from the free neutron 
decay process). If |V| dominates over |V|, neutrino amplitude evolution becomes non- 
adiabatic. 

The ratios of the first three rates are: 



v, 



H 



r ( t 



(ii) 



H VMeV, 

In these ratios, an average E = 3.15T is assumed. Since we are only concerned with 
|5m 2 | > 0.1 eV 2 , as long as sin 29 > 10 -6 (the minimal mixing required to amplify L Ua , see 
Eq. (0)), the neutrino transformation rate easily dominates over the Hubble expansion and 
the weak interaction at T ~ 1 MeV. 

We always have V x = HV X . At T < T rcs /2 and away from resonances so that we can 
assume Vj < \5m 2 \/2E ~ V£ ~ V z , we have \V Z \ ~ \(H - L u jL Va )V z \ ~ \HV Z \. (Note 
that L Va is limited to the T~ 4 growth in Eq. (p.) Therefore at T < T rcs /2, the active- 
sterile neutrino transformation channels in our problem can be treated as ordinary oscillation 
channels with adiabatically varying mixing parameters except possibly at the resonances. 
We will discuss the question of adiabaticity at resonances later. 

At T < T res /2 when we can neglect V<[ , the resonance condition V z = — 5m 2 /2E cos 2$ + 
V£ = gives the asymptotic values of L Va in Eq. (0). Note that for L Ua < [L Va > 0) the 
v a [pa) transformation channel is matter-enhanced. The fraction F of resonantly converted 
v a ipa) in the total v a (P a ) distribution is 

6 



V, 



de 



dV z 



2e res sin 26>/(e res ) 



(12) 



where e = E/T, /(e) « 2e 2 /3£(3)(l + e e ) and e res satisfies I4( e res) = 0. The resonance e 
remains stationary (or varies very slowly) and L Ua oc T -4 if 

3„ 3 



-e res sin 26f(e TCS ) ~ 0.1 sin 20. 



(13) 



In general, the stationary e res has a value ~ 0(0.1) @,|]- When \L Va \ > 0.1 sin 26 1 , |L„J 
cannot keep up with the T~ 4 growth. Consequently e res must increase as e res oc T~ A \L Va \~ l 
at low temperatures when V T can be safely neglected. 

The resonance is adiabatic when the resonance region is characterized by slowly changing 
e res . As the resonance sweeps through the neutrino energy spectrum, a complete conversion 
of v a to u s at the resonance is possible only if the following adiabatic conditions are met: 



V 2 

y X 




V x sm29' 


VI 

v z 




H 



10 y 



1 5m 2 

~W 2 



T 



MeV 



sin 2 26 > 1, 



(14) 



and 
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de 



dV z 



de r 



dL, 



>V X 



de 



dV z 



(15) 



The first condition simply implies that the timescale of completing the resonance has to 
be much longer than the neutrino oscillation period at resonance. This is satisfied if 
\5m 2 \sm 2 29 > 10~ 9 eV 2 . The second condition requires that e res move slowly through 
the spectrum. This amounts to 



Li/ a ^ ^^res/(^res); 



(16) 



which is satisfied for e res 's that cover the bulk of the neutrino energy spectrum. Adiabaticity 
at resonances is therefore a valid assumption for u a <-> v s mixing with vacuum mixing angles 
which are not too small but which are still within the BBN bound Eq. ([|). 

Given adiabaticity, the growth of L Ua (assumed to be positive for simplicity and no loss 
of generality) as e res sweeps through the neutrino energy spectrum can be easily estimated 
by solving 



3 r e ^(T) 

L Vot (T)^-J o Pf(e)de (17) 

and the resonance condition V^(e res ) = which at low temperatures is equivalent to 

^ |<5m 2 /eV 2 | 
ercsUJ- i6(T/MeV) 4 L, Q ' 1 J 

In Eq. ( |T7|) /3 takes account of the effect of collisions that redistribute energy among neutri- 
nos. 

We can identify two extreme cases. When the collisions are too inefficient to change 
the neutrino distribution at e > e res (such as when T < 1 MeV), j3 = 1. In another limit, 
8L Ua /3 which obtains when the collisions are highly efficient (such as when T > 1 
MeV) and neutrinos are always distributed thermally. Eqs. (|T7[) and ([T8|) give a solution 
in fair agreement with the results obtained by solving Eq. (18) of Foot and Volkas ||. In 
Figure 1 we plot our results in terms of L Ua vs. m 2 , — m 2 s at various temperatures. From 
Figure 1 we can deduce a power law relation applicable to \L Va \ < 0.1, 

\L Va \ « 0.05|5m 2 /eV 2 | 2/3 (T/MeV)~ 8/3 . (19) 

We note that this power law applies only in the stage when the resonance sweeps through 
the neutrino energy spectrum. When the resonance is stationary (Eq. fljjD ), the dependence 
is \L Va \ oc \8m 2 \T~ A instead. 

The asymmetry L Va is generated as the resonance conversion region moves up through 
the neutrino energy spectrum. This suggests a distortion of the v a (or v a ) energy spectrum 
(see also Ref. jnj). Indeed, when L Va < 0.1, most of the L Ua is generated at the lowest 
temperatures (Eq. []T9[ ) , when the neutrino scattering processes that tend to thermalize the 
neutrino spectrum are the most inefficient. The fact that the resonant transformation of v a 
to v s (or v a — > v s ) starts at lower neutrino energies only further deepens the inefficiency of 
neutrino re-thermalization, as neutrino interaction cross sections scale roughly linearly with 
neutrino energies. In Figure 2, we plot a semi- analytical calculation of the v a spectrum at 
T = 1 MeV for a v a — > v s resonant transformation (which generates a negative L Ua ) with 
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ml a — m 2 Vs = 1 eV 2 . In the calculation, the thermalization process is approximated as a 
relaxation process (with a rate T) driving the system toward a thermal distribution. As a 
result of the inefficiency of this process, for cases with \Sm 2 \ < 1 eV 2 (so that L Va < 0.1 at 
T ~ 1 MeV), the u a neutrino spectrum at and below its thermal decoupling temperature 
T ~ 1 MeV can be well approximated by a thermal spectrum with a low energy cut-off. The 
u a deficit below the cut-off energy results in the L Ua asymmetry. In the mean time, u a is 
not subject to resonant transformation. Its spectrum is therefore not significantly changed, 
due to the inefficiency of neutrino pair production. (The opposite is true if L Ua > 0: the u a 
distribution will have its lower energy region truncated but the v a distribution will remain 
intact.) 

III. DIRECT GENERATION OF ELECTRON-NEUTRINO ASYMMETRY BY 
RESONANT v E v s TRANSFORMATIONS 

If a = e (the direct scheme to generate an asymmetry in the z/ e z/ e sector), the v e or V e 
spectral distortion will directly impact the neutron-to-proton ratio at the weak freeze-out 
temperature, and hence the 4 He yield. Figure 3 shows the changes in n<->p rates due to the 
neutrino spectral distortion in the case m 2 e — m 2 s = 1 eV 2 . When L„ e > (i.e., a deficit of 
low energy P e ), the major effect is an enhanced neutron decay rate at low temperatures due 
to the reduced Pauli-blocking of P e . For reaction p+z/ e ^n+e + , the low energy deficit in V e 
is of little significance because only v e with E > 1.9 MeV can participate in the reaction. 
Conversely, its reverse reaction mostly generates P e at the higher end of the energy spectrum. 
Its rate is therefore insensitive to the spectral distortion at the lower end. When L Ue > (a 
deficit of low energy v e ), the rate for n+z/ e ^p+e~ is significantly reduced while the reverse 
rate is slightly increased. Figure 4 shows the resultant AY from the spectral distortion as a 
function of m 2 e — m 2 s for both L Ue > and L Ue < 0. The disparity between the two cases 
of opposite L Ue is transparent from Figure 3: the change in n^p rates is much larger when 
L Ve < 0. 
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The v e Majorana mass limit (which is uncertain by a factor of 2, ranging from m Ve < 0.45 
|T7f ) implies an upper limit m 2 ^ —vr? v < 1 eV 2 . [] Figure 4 shows that 
the maximally allowed reduction in Y is only « 0.0021, about 1% of the standard prediction. 
But the maximally allowed increase of Y can be as high as ~ 0.022, a 9% effect. An increase 
this large in the predicted primordial 4 He abundance would have already been too large to 
accommodate observations |7|||. 



eV to m Ve < 1 eV [15 



IV. INDIRECT GENERATION OF ELECTRON-NEUTRINO ASYMMETRY BY 

NEUTRINO TRANSFORMATIONS 

For convenience in the indirect scheme, we assume that a L Ut is first generated by a 
v T <-> v s transformation process. This asymmetry may then be transferred to L Ue by a 
resonant v T <-» v e oscillation ||. (Ordinary oscillations without resonance cannot transfer the 
asymmetry efficiently) In fact we are likely to have Sm 2 TS ^ = m 2 Us —m 2 UT « 5m, Te \ = vn 2 Ve —w? Vr 
if m„ T ^> m Ue , m Us (which will be the case in order to have an appreciable impact on the 
primordial 4 He abundance). Such a neutrino mass spectrum is consistent with a simultaneous 



solution of the solar neutrino problem [ITSJ and the atmospheric neutrino problem [I9]|. For 



the moment, we will simply assume 5m 2 TS ^ ~ 5m 2 re y 

The matter asymmetry contributions to the effective potentials of the two neutrino trans- 
formation channels become 

Vfc s) w ±0.35G F T 3 (2L l/r - L Ve ), Vfa w ±0.35G f T 3 (L 1/t - L Ve ). (20) 

Apparently, the v T <-> v s and ^ T «-> z/ e resonances when > (or the z/ T <-> f s and z/ T <-> v e 
resonances when L Ut < 0) do not simultaneously share the same part of the neutrino energy 
spectrum. 



1 We point out that it might be possible to construct mixings between left-handed electron neu- 
trinos and left-handed sterile neutrinos. In such a case the Majorana mass limit will not apply but 
the standard weak interaction theory will then need to be reconsidered. 
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Guaranteed adiabaticity (i.e., satisfying Eq. fl4j), the efficiency of resonant neutrino 
conversion is still determined by whether or not the neutrino collision time scale dominates 
over the timescale for a complete v T to v e (or v T to u e ) transition at resonance (the resonance 
width). The collision timescale is important because the two resonances in Eq. ( pop do not 
overlap, so for example, any deficit in v T caused by the v T <-> v s resonant transition relies 
on neutrino scattering to redistribute neutrinos into the energy region where the v T <-> v e 
resonance occurs. In previous work [f| this redistribution has been assumed to be instant. 
However, this is not a good approximation at T < 5 MeV as we will show below. 

The ratio of the resonance timescale to the collision timescale is 



r ~ ]-sin2# (Te) ~ (^) 3 sin2^ (Te) . (21) 



Calculations based on active-active neutrino transformation in Type II supernova nucleosyn- 
thesis suggest that sin 2 29 ( Te ) < 10~ 4 for m 2 r — m 2 e > 1 eV 2 PCfl , and the Bugey experiment 
constrains sin 2 26 l ( Te ) < 0.04 for m 2 r — m 2 e < 1 eV 2 |[21|| . Therefore, the collision timescale 
dictates the growth of L Ve at T < 5 MeV for \6m^ TS ^\ rs |5m 2 re) | > 1 eV 2 . 

Above 5 MeV, the collisions may be deemed instantaneous, and the equations of ref. 
H become valid. But a side effect of generating a large L Vt above ~ 5 MeV by the v T <-> 
v s mixing (by having m 2 r — m 2 s > 10 4 eV 2 ) is bringing sterile neutrinos into chemical 
equilibrium ||. As a result, even though a reduction as large as AY = —0.006 may result 
from a positive L Ve alone in the indirect scheme |J, the extra sterile neutrinos produced will 
increase Y by at least as much by bringing v s into chemical equilibrium through the v T <-> v s 
resonant transformation. The net effect is an increase of Y instead of a reduction, even if 
L Ve is positive. 

We again model the transfer of L Vt to L Ue by assuming that active neutrinos tend to 
their equilibrium distributions with a rate T. Figure 5 shows the growth of L Vt , L Ue and 
the accompanying increase in the total neutrino energy density (including that of sterile 
neutrinos) in the units of that of one unperturbed active neutrino flavor, for one particular 
set of mixing parameters. Figure 5 shows that L Ur indeed grows initially according to 
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Eq. (0), and tapers off at L Vt pa 0.22 when most of v T or v T have undergone resonances. 
The transfer of L Ut to L Ue is efficient at T > 2 MeV, but freezes out below ~ 2 MeV. It 
freezes out at a higher temperature than L Vt because the v T -v e (or v T -v^) resonance occurs 
at a higher energy than the v T -v s (or v T -T>s) resonance (Eq. ||20|| ). The resonance therefore 
sweeps through the v e (or i> e ) spectrum faster. The figure also shows that increase in the 
total neutrino energy density in this case (about 2%, or AN U ~ 0.07) is moderate. 

The spectra of u e and v e in the indirect scheme is only slightly distorted. Figure 6 shows 
the modified v e spectrum when L Ut , L Vr > for m 2 T — m 2 e = m 2 r —m 2 Us = 100 eV 2 , compared 
to an unperturbed active neutrino spectrum with zero chemical potential. The distortion is 
small because the transfer of L Ut into L Ue occurs in the entire energy distribution (albeit at 
different temperatures), unlike in the direct scheme when the generation of L Ue occurs only 
in low energies. The distortion in v e and v e spectra can be well approximated by an overall 
multiplication factor 1 + 5±. The net asymmetry is therefore L Ue = 3(5 + — 5-)/8, and the 
percentage increase in the v e v e number density due to pair production is (5 + + 5-)/2. 

By modifying the standard BBN code with the new v e and v e spectra, and with the 
increased total neutrino energy density, we obtain their effects on Y in Figure 7. At m 2 ^ — 
ml, = mf, — m? Us 100 eV 2 the effect on Y is dominated by the asymmetry in the v e v e 
sector. But as m 2 Vr — m 2 e = m 2 T — m 2 g increases, the increased total neutrino energy 
density gradually becomes significant. The increased energy density causes Y to increase, 
regardless of the sign of the neutrino asymmetry. As a result of these two factors, a maximal 
reduction AY pa —0.005 is achieved in cases of positive lepton number asymmetries when 
m 2 T — m 2 e = m 2 T — m 2 s ~ 100 to 300 eV 2 . This mass-squared-difference, however, implies 
that tau neutrinos are unstable, based on cosmological structure formation considerations 

Figure 7 is very different from the previous estimates of Foot and Volkas . For example, 
Foot and Volkas have argued for a possible reduction AY pa —0.006 across the mixing 
parameter range 10 < m 2 T — m 2 e = m 2 T — m 2 s < 3000 eV 2 . While our result clearly shows 
a concave feature of AY in this range, with a maximum at pa —0.005. Foot and Volkas' 
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result also indicated that AY is smaller in the positive direction (when L Ur , L Ue < 0) than 
in the negative direction (when L Ut , L Ue > 0). Our analysis indicates the opposite: when 
L„ T , L Ue < 0, the changes in Y due to the spectral asymmetry and the extra neutrino energy 
add constructively; while L Ut) L Ve < 0, these two effects add destructively. AY, therefore, 
is larger in the positive direction than in the negative direction. 

These differences, we believe, stem from our detailed consideration of neutrino spectrum 
distortion and its time dependence. These factors are crucial to the neutron-to-proton 
freeze-out process, and in turn the primordial 4 He yield. 

V. SUMMARY 

We have calculated the spectral distortions for neutrinos and the time dependence of 
the neutrino distribution function during the lepton asymmetry generation via resonant 
active-sterile neutrino transformation. We have included these crucial effects in our BBN 
calculation assessing the effect on the primordial 4 He abundance of the possible lepton num- 
ber asymmetry in the v e v e sector. We conclude that the possible increase in the primordial 
4 He yield, as a result of a negative lepton number asymmetry, can be substantial. The 
maximal increase can be as high as ~ 0.01 to 0.02 (or 5 to 9%) for mixing parameters 
that are consistent with neutrino mass constraints. The possible decrease due to a positive 
lepton number asymmetry, however, is limited to < 0.002 (or < 1%) if the asymmetry is 
generated by a resonant v e <-> v s mixing, or < 0.005 (or < 2%) if the asymmetry is gener- 
ated by a three-family resonant mixing scheme. The magnitude of these possible changes 
in the primordial 4 He abundance induced by the neutrino-mixing-generated lepton num- 
ber asymmetry is comparable to or greater than the uncertainty of current primordial 4 He 
measurements. Therefore, the role of resonant active-sterile neutrino mixing in Big Bang 
Nucleosynthesis cannot be underestimated. 

X. S., G. M. F. and K. A. are partially supported by NSF grant PHY98-00980 at UCSD. 
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Figure Captions: 

Figure 1. The magnitude of lepton asymmetry as a function of 5m 2 at various tempera- 
tures. The bands denote the range of the asymmetry enclosed by the two extreme cases: 
(1) collisions are completely inefficient (upper limits); (2) collisions are completely efficient 
(lower limits). 

Figure 2. The solid curve: the calculated v a (a = e, //, r) distribution function. The dashed 
curve: an unperturbed thermal neutrino distribution function with zero chemical potential. 

Figure 3. The change in n^p rates due to u e (if L Ue > 0) or v e (if L Ue < 0) spectral distortion 
for m?. —ml =1 eV 2 . 

Figure 4. The impact on the primordial 4 He abundance Y if an asymmetry in u e u e is 
generated by a resonant v e <-> v s mixing in BBN. Baryon number density to photon number 
density ratio is set to rj = 5.1 x 1CT 10 . 

Figure 5. The solid curves: \L Vt \ and \L Ue \ as a function of temperature for m 2 T — = 
ml T — m 2 Vs = 100 eV 2 . The dashed curve: the increase in the total neutrino energy density 
as a function of temperature, normalized by the energy density of one thermalized active 
neutrino flavor with zero chemical potential. 

Figure 6. The solid curve: the calculated v e distribution function in the indirect mixing 
scheme. The dashed curve: an unperturbed thermal neutrino distribution function with 
zero chemical potential. Inset: the ratio of the two distribution functions vs. neutrino 
energy. 

Figure 7. The impact on the primordial 4 He abundance Y in the indirect neutrino mixing 
scheme, as a function of m 2 T —m 2 Ue = —nT>l 3 - Baryon number density to photon number 
density ratio is set to r\ = 5.1 x 10~ 10 . 
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